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ABSTRACT 

V838 Mon erupted at the beginning of 2002 becoming an extremely luminous star with L ~ 10^ Lq. 
Among various scenarios proposed to explain the nature of the outburst the most promising is a 
stellar merger event. In this paper we investigate the observational properties of the star and its 
surroundings in the post outburst phase. We have obtained a high resolution optical spectrum of 
V838 Mon in October 2005 using the Keck I telescope. We have identified numerous atomic features 
and molecular bands present in the spectrum and provided an atlas of those features. In order to 
improve the spectrum interpretation we have performed simple modeling of the molecular bands. 
Our analysis indicates that the spectrum is dominated by molecular absorption features arising in 
photospheric regions with temperatures of ~2400 K and in colder outer layers, where the temperature 
decreases to ^500 K. A number of resonance lines of neutral alkali metals are observed to show P-Cyg 
profiles. Particularly interesting are numerous prominent emission lines of [Fe II]. All of them show 
practically the same profile, which can be well described by a Lorentzian profile. In the blue part of 
the spectrum photospheric signatures of the B-type companion are easily seen. We have fitted the 
observed spectrum with a synthetic one and the obtained parameters are consistent with the B3V 
type. We have also estimated radial and rotational velocities of the companion. 

Subject headings: binaries: spectroscopic — stars: emission-line — stars: individual (V838 Mon) — 
stars: mass loss — stars: winds, outflows — stars: variables: other 



1. INTRODUCTION 

The eruption of V838 Mon was discovered in the be- 
ginning of January 2002. Initially thought to be a nova, 
the object appeared unusual and enigmatic in its nature. 
The eruption, as obse rved in the optical, lasted about 
three months ([Mu nari et al.l l2002bl : 1i<imeswenger et all 
12002'; Cra use et al.l 2003). After developing an A-F su- 
pergiant spectrum at the maximum at the beginning of 
February 2002, the object evolved to lower effective tem- 
peratures and in April 2002 it practically disappeared 
from the optical, remaining very bright in the infrared. 
At the same time a B3V companion to the erupted object 
was discovered in the optical ( Munari et al. 2002a). A 
detailed analysis of the evolution of the object in the ou t- 
burst and decline can be found, e.g. in lTvlendal (f2005f) . 

Several mechanisms have been proposed to explain 
the eruption of ¥ 838 Mon, including an unusu al nova 
(llben fc Tutukovl [19921. a lat e He-shell flash (iLawloH 
I2005D, and a stellar merger (Soker & Tvlenda '200%^ 
They have critically been discussed in iTylenda & Soken 
(|2006). These authors conclude that the only mechanism 
that can satisfactorily account for the observational data 
is a collision and merger of a low-mass pre-main-sequence 
star with a ~ 8 Mq main-sequence star. 

In October/November 2006 the B3V companion sig- 
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nificantly faded and a strong Ha emission appeared 
in th e spectrum of V838 Mon (jCoranskiil l2006t iBondl 
l2006f l. Late in 2004 an emission-line spectrum, com- 
posed mainly of [Fe II] lines, started to develop 
(jBarsukova et aP l2006n , re aching its maximum around 
the 2006 eclipse-like event (jMunari et al.ll2007[ ). 

In the present paper we present and discuss a high res- 
olution spectrum of V838 Mon acquired with the Keck I 
telescope in October 2005. At that time the [Fe II] 
emission-line spectrum was already well developed and 
the object was a year before the eclipse of the B3V com- 
panion. 

2. OBSERVATIONS AND DATA REDUCTION 

A high resolution spectrum of V838 Mon was obtained 
on 2005 October 13 UT with the High R esolution Echelle 
Spectrometer (HIRES. IVogt et al.|[l99l attached to the 
Keck I telescope. We used the C4 decker, which pro- 
vides a slit with dimensions 1'.'148 x 3'.'5. The resultant 
resolving power was of i? ~ 34 000. The HIRES instru- 
ment makes use of a 3-chip mosaic as a detector. We 
obtained two exposures of the object with a total time 
of 960 s. 

Data reduction was performed with the IRAF'' stan- 
dard procedures (Masscv 1997). Data were debiased and 
flat-fielded using the ccdproc task. Cosmic ray events 
were remo ved from the C CD fr ames with a method de- 
scribed in Ivan DokkumI (|200l[ ). After a careful check 
of the removal we are confident that the procedure did 
not introduce any spurious effects to the spectra. Or- 
der tracing and extraction were done using the tasks in 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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the echelle package. The orders corresponding to the 
bluest part of the spectra were underexposed and we 
were not able to extract them, even with a pinhole trace 
as a reference. The wavelength calibration was based 
on ThAr lamp observations and it is accurate to within 
0.003 A. The two extracted echelle spectra were aver- 
aged giving a spectrum with a total wavelength coverage 
3720 - 7962 A. In the range there are, however, numer- 
ous gaps caused by the inter-chip spaces on the mosaic 
and the orders that do not overlap. 

Using an extracted pinhole frame we corrected all the 
orders of the echelle spectrum for a blaze effect. Inspec- 
tion of the overlapping and blaze-corrected orders re- 
vealed that, although correctly flattened, they were sys- 
tematically tilted in such a way that the long- wavelength 
end of an order was always higher on a relative intensity 
scale than the short-wavelength end of the neighboring 
order. The discrepancy was of 15% of the average in- 
tensity of the order. We reduced the tilt by multiplying 
each order by a linear function (i.e. with 1 in the middle 
of an order range, and 1 ± 0.075 at the edges). All the 
overlapping orders were then safely merged and, subse- 
quently, we got a spectrum composed of three parts: (i) 
completely covered wavelength range over 3720 -4883 A, 
(ii) wide range over 4951 — 6469 A with two narrow gaps 
at 6308.5 - 6310 A and 6401 - 6424.5 A, and (lii) a part 
of the spectrum over the range 6544 — 7962 A with 10 
non-overlapping orders. The three pieces we call here- 
after parts blue, green and red, respectively. 

Since no spectrophotometric standard stars were ob- 
served during our observing run at Keck I, we decided 
to flux calibrate the HIRES spectrum o n a base of the 
method described in lSuzuki et al.l (|2001h . namely by us- 
ing a flux calibrated spectrum of V838 Mon obtained 
elsewhere with a lower spectral resolution. Such a low 
resolution spectrum was acquired on 2005 November 22 
with the Grating Spectrograph with SITc mounted at the 
1.9 m Radcliffe telescope at the South African Astronom- 
ical Observatory (SAAO). We made use of a grating #7 
with 300 lines mm~^, and a slit with a projected width 
1'.'5. The last figure was comparable to seeing conditions 
during the observation (~1'.'4). The total duration of the 
exposure was 1000 s. Three spectrophotometric standard 
stars, i.e. LTT9239, EG21, and LTT3218, were observed 
during the night. This allowed us to calibrate the low res- 
olution spectrum in absolute flux units, while the wave- 
length calibration was performed with ThAr reference 
lamp spectra. All the data reduction and calibrations 
were carried out with standard IRAF procedures. The 
extracted and flux calibrated spectrum covers the range 
3817 - 7235 A with a resolving power of i? ~ 1 000. 

The flux calibration of the HIRES spectrum was per- 
formed as follows. The high and low resolution spectra 
were corrected to the heliocentric rest frame and they 
were cross-correlated in order to find a relative shift be- 
tween them in the wavelength domain. We found a shift 
of 5.5 A (note that the value is close to the spectral reso- 
lution in the SAAO observations), which was eliminated 
by displacing the low resolution spectrum. Next, the 
HIRES spectrum was smoothed to the resolution of the 
observations at SAAO and the spectra were divided one 
by the other. Rather than to calibrate each single order 



of the echelle spectrum, as it was done in iSuzuki et al.l 
((2001I) . we decided to calibrate separately only the three 
mentioned above parts of the merged spectrum. To each 
of the pieces a low order polynomial was fit ted giving a 
conversion ratio (CR, cf. ISuzuki et a l. 2001). Since the 
SAAO spectrum has a narrower wavelength range than 
the high resolution one, we extrapolated the CRs to the 
full ranges of the blue and red parts. Finally, the HIRES 
spectrum was put on the absolute flux scale by dividing 
its parts by the appropriate CRs. We note that the flux 
calibration of our final spectrum is uncertain, particu- 
larly in the ranges not covered by the SAAO observations 
and in the regions where strong atmospheric molecular 
bands are present (the object was observed at unequal 
air masses in the different facilities). 

Although the low resolution spectrum was obtained 
about six weeks after the HIRES observations, we are 
quite confident that the spectrum of V838 Mon did not 
change significantly during this period in the observed 
spectral range. A comparison of the spectra smoothed 
to the same resolution shows that there are no signifi- 
cant discrepancies which would affect the flux calibra- 
tion. Moreover, from the photometric behavior of the 
object, illustrated e.g. on the V. Goranskij's web page^, 
one can see that the assumption of a constant flux in B 
and V bands during the period is justified. 

The fiux calibrated spectrum was then corrected for 
the interstellar redd ening using t he IRAF's task deredden 
with Eb-v = 0.90 (Tylenda 2005) and with a standard 
ratio of total to selective extinction Ry = 3.1. 

3. THE SPECTRUM 

The fiux-calibrated and dereddened spectrum of 
V838 Mon is presented in Figs. [THSl For a comparison a 
synthetic spectrum of a B-type star (see Sect. E]) is over- 
ploted with the blue line. Most of identified atomic and 
molecular features are indicated in the figures. Wave- 
lengths are given in the heliocentric rest frame and the 
units of the flux are 10~^^ erg s~^ cm~^ A~^. 

Our observations reveal an unusual and very complex 
spectrum of V838 Mon. It contains contributions from 
several physically different environments. In its blue part 
a photospheric spectrum of the B3V companion can be 
clearly seen, while the green and red parts are dominated 
by molecular absorptions that severly affect the flux of 
a very cold supergiant. A particularly striking feature 
of the presented spectrum is the presence of prominent 
emissions of [Fe II]. Finally numerous, mostly resonance, 
lines showing P-Cyg profiles as well as a few pure ab- 
sorptions of neutral atoms can be easily found. 

In the case of such a complex spectrum, detailed iden- 
tification of all the features is a difficult task. The identi- 
fication procedures and their results are described in the 
following sections. We also present results of basic mea- 
surements, e.g. radial velocities of different line systems 
and molecular bands. 

4. ATOMIC LINES 

The atomic line identification was based mainly on the 
Atomic Spectra Database Lines Form^ (provided by the 
National Institute of Standards and Technology, NIST), 

^ http://jet.sao. ru/~goray/v838mon . htm 
http: //physics .nist . g ov/PhysRef Data/ AS D/llnes_f orm . html | 
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Fig. 4. — Profile of the [Fe II] lines. Black solid line: the profile 
averaged out from 7 lines (see the text for more details). Red 
dotted line: the best least-squares fit of a Lorentzian profile to the 
average profile (FWHM of 83.6 kms-\ center at 13.0 kms~i). 
Blue dashed line: local continuum. 

the Atomic Line List v ersion 2 . 04 by van Hoof^, and on 
the multiplet tables bv lMoord (| 19451 ). 

All the identified atomic emission and absorption fea- 
tures are listed in Tables[l]and[5] (respectively) in order of 
their laboratory wavelengths. Columns (l)-(3) of the ta- 
bles contain: labor atory wavelen gth, ion name, and mul- 
tiplet number from lMoori (|1945l ). Column (4) in Table [1] 
gives integrated (dereddened) fluxes of the emission lines 
(for which a reliable measurement was possible). In the 
same column of Table[2]the profile type (pure absorption 
or P-Cyg) is specified. Notes on the features are given in 
the last column of the tables. 

4.1. Emission lines 
4.1.1. Forbidden lines 

A great majority of the emission features seen in the 
spectrum are forbidden lines of Fe II. They are remark- 
ably numerous. We have identified 49 lines from 11 differ- 
ent multiplets. All of them arise from levels with energies 
(above the ground level) between 2.0 and 3.4 eV. Most of 
them are very strong and often blended with other spec- 
tral features. In addition to these forbidden transitions, 
we have also found two permitted lines of Fe II at 5018 A 
and 5169 A (a third line from the same multiplet falls in 
a gap between two chips). 

All the unblended [Fe II] emissions show practically 
the same but rather unusual shape. It can be well fit- 
ted with a theoretical Lorentzian profile, as can be seen 
from Fig. (4] The profile shown in this figure is a mean 
shape obtained from seven strongest and "clean" [Fe II] 
lines. When averaging, the normalized profiles of indi- 
vidual lines were weighted according to the rms^ values 
of the local noise in the spectrum. From 30 [Fe II] lines, 
for which we successfully fitted a Lorentzian profile using 
IRAF's splot, we have obtained a typical full width at half 
maximum^ (FWHM) of 76.3±17.7 km s"^ (median and 
standard deviation). All the lines have practically the 
same radial velocity. Measurements performed with the 

Ihttp : //www ■ pa . uky . edu/~peter/ato mlc/ 1 

* All line widths given in this paper are not corrected for the in- 
strumental profile, which has an intrinsic width of about 9 km s~^. 



IRAF's rvidlines task gave a heliocentric radial velocity 
of T4=13.3±0.7 km s~^ (mean and standard deviation). 

Other forbidden lines seen in our spectrum include 
those of [O I] 5577 A, 6300 A, 6364 A and [N I] 
5198 A. These are rather weak features, observed in spec- 
tral regions strongly contaminated with molecular bands. 
Their measurements are therefore uncertain. The mea- 
sured FWHM and radial velocity of the [O I] lines are 
80 ± 30 km s~i and 40 ± 10 km s"\ respectively. The 
[N I] line has a radial velocity of about 80 km s~ . 

It should be noted that some of the emission lines listed 
in Table [1] have uncertain identification. Two strong fea- 
tures at about 5041 A and 5056 A were identified as Si II, 
but other identifications are possible. Three lines recog- 
nized as [Ni II] 3993 A, 4201 A, and 4326 A, are weak, 
which makes their identification uncertain. However, the 
ionization potential of Si and Ni being similar to that of 
Fe and low excitation energies of the proposed multiplets 
make the identification probable. 

4.1.2. Balmer lines 

Inside the cores of the photospheric Balmer lines of the 
B-type companion weak emission features can be found. 
Their presence is evident when one compares the ob- 
served photospheric profiles with a synthetic spectrum 
of an early B-type star (see Fig. [1]). The following dis- 
cussion is limited to H/3 only, as the emission feature is 
most evident in this line. In order to better characterize 
the emission we have subtracted the underlying stellar 
absorption. The procedure was as follows. To the wings 
of the absorption line we fitted a profile in the form: 

0(A) =dexp{[-a(A±Ac)'' + c]-i} (1) 

with a, 6, c and d as free parameters. The central wave- 
length, Ac, was taken as the laboratory wavelength of 
the H/3 line corrected for the radial velocity of the B3V 
star (see SectJB]). Eq. ([1]) is based on an empirical pro- 
file given in IChauville et all (|2001f) . However, for our 
purpose it was generalized to treat also the local con- 
tinuum level as a free parameter. The fitting procedure 
used the least-squares method and it was applied to all 
the data poi nts in the H/3 w i ng [i.e ., not only to three 
points, as in IChauville et ahl (|20011]. The core region 
dominated by emissions was, obviously, excluded from 
the fit. The extracted emission feature is shown in Fig.[5l 
It can be interpreted as a single emission with an absorp- 
tion superimposed on it. In this case, fitting Lorentzian 
profiles gives the emission component with a FWHM of 
~42 km s~^ at a radial velocity of ~14 km s""'^, while 
the absorption component appears at ~10 km s~^ and 
has a FWHM of ^^24 km s~^. Note that the parameters 
of the emission component in this case are close to those 
of the [Fe II] lines. Alternatively, if the observed profile 
is fitted with two emission components, they have radial 

velocities of 10 and ~28.8 km s"\ and FWHMs of 

~30 and ~24 km s~^, respectively. 

4.2. P-Cyg features 

The spectrum displays a number of strong resonance 
lines showing P-Cyg profiles. The most prominent ones 
are those of Mn I 5395 A and 5433 A, Cr I 5410 A, Ba I 
5535 A, Ca I 6573 A, K I 7699 A, Rb I 7800 A (only 
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Fig. 5. — Emission feature extracted from the H/3 photospheric 
absorption line of the B3V companion. The spectrum was signifi- 
cantly smoothed from the original resolution via boxcar with a box 
size of 15 pixel. 

the Cr I line is not a resonant one). They are shown in 
Fig. [HI Other hnes showing P-Cyg signatures, although 
not as notable as in those listed above, are Mg I 4571 A, 
Fe I 5110 A, and possibly also a few weak lines of Cr I. 

Most of the P-Cyg profiles have a deep absorption com- 
ponent extending from ~50 km (as a cross-over from 
emission to absorption) to ^ —95 km s~^. The most ex- 
treme exception is the K I line, which shows blueshifted 
absorption down to —145 km s^^. Emission compo- 
nents peak at ~80 km s^^. 

The shapes of the absorption components indicate in- 
homogeneities in the outflowing matter. Particularly in 
the profiles of Rb I and Mn 1 one can clearly see an ad- 
ditional narrow absorption component (NAC) centered 
at —82 km s~^ [see panel (b) in Fig. [B]. The measured 
FWHM of the narrow feature inside the Rb I 7800 A line 
is 19.3 km s^^. 

As noted above, the Mg I 4571 A line shows only weak 
signatures of a P-Cyg profile. Indeed, it is dominated by 
a slightly asymmetric emission and exhibits only a very 
weak absorption at approximately the same velocity as 
the NAC. 

The profile of Fe I 5110 A is also exceptional. Its ab- 
sorption component is so deep that it removes all the 
light from the cold star and reaches the flux level of 
the B3V companion (see top panel in Fig. [2]). In addi- 
tion, it is considerably broader than most of the other 
P-Cyg absorptions and extends from ~80 km s~^ to 
~ —120 km s~^. Most probably the A5110 absorption 
is saturated. Other Fe I lines identified in the spectrum 
appear as absorptions only, without any signs of emission 
components (see Sect. 14. 3p . 

As Table [2] indicates, almost all the abundant alkali 
metals show P-Cyg profiles. Thus, one would expect 
that it should also be the case of the Li I doublet at 
6708 A. However, in the raw spectrum no clear emission 
component can be seen close to the clear absorption of 




1.8 
1.6 

1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 

1.6 
1.4 

1.2 
1 

0.8 
0.6 
0.4 
0.2 


1.2 



0.8 



1 - 



I I I I I I I I I ) I I I I I I I I I I I 




.(c) , 

I I I I I I I I I I I M 




-200 



-100 







100 



200 



300 



Heliocentric velocity (km s ) 



Fig. 6. — Sample of atomic spectral features. Before the regions 
were individually normalized to the local continuum, the flux from 
the B3V star had been subtracted, (a) Examples of the most 
prominent P-Cyg lines. The spectra of Cr I, Ba I and Fe I, were 
smoothed via boxcar 9. The K I profile is contaminated by the ISM 
absorptions, (b) P-Cyg profiles with a narrow absorption compo- 
nent (NAC). The normalized fluxes of the Mg I were divided by 
1.7 and the spectrum was smoothed via boxcar 13. The Mn I lines 
were smoothed with boxcar 9. (c) Double absorption components 
of three Ti I transitions; the upper and lower plots were vertically 
shifted by ±0.55. (d) Two absorption lines of Ca I. The vertical 
lines along the panels mark the velocity of the NAC (dashed) and 
of the source of the outflow (dotted) as estimated from molecular 
bands (see Sect. \5A\ . 
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Li I. The doublet falls in a region that is dominated by a 
strong absorption band of TiO and it is difficult to define 
an underlying molecular baseline. Moreover, the doublet 
resides very close to the 7 (1,0) bandhead of TiO and 
most probably a part of the emission component of Li I 
is absorbed by the saturated molecular band. A similar 
situation also occurs for the Rb I 7948 A line. Its profile 
is strongly affected by the (B— X) (0,0) absorption band 
of VO and the expected emission component cannot be 
seen on the raw spectrogram. 

4.3. Absorption lines 

Another group of identified lines are pure absorptions. 
Beside the interstellar features (see Sect. [7|) and photo- 
spheric lines of a B-type star (see Sect. [S]), absorptions 
of Fe I 5060 A, 5128A, 5247 A, 6551 A, Ca I 6102 A 
and 6122 A, and numerous Cr I lines were found. Some 
of these lines are in the same multiplet with a line that 
exhibits a prominent P-Cyg profile. Thus, it is quite 
possible that in a more sensitive observations the lines, 
classified here as absorptions, would appear as P-Cyg 
profiles. Among the absorptions the strongest is the Fe I 
5060 A line, which, similarly to the P-Cyg absorption 
component of Fe I at 5110 A (see Sect. 14.21) . appears as 
a very broad feature and, probably, is saturated. All the 
absorptions, except the A5060 line, are narrower than the 
absorption components in the P-Cyg profiles and always 
fall in the velocity range defined by P-Cyg absorptions, 
indicating that they originate in the same parts of the 
outflowing gas. 

The absorption lines of Ti I found at 5147 A, 5152 A, 
and 5426 A are of special interest. They arise from lev- 
els with energies of 2.3—2.4 eV. In our spectrum the hues 
consist of (at least) two well separated components, what 
is clearly seen in panel (c) of Fig. \6\ The redshifted ab- 
sorption component appears at ~75 km s"^, close to the 
peaks in the P-Cyg emissions. The blueshifted sharp fea- 
ture, seen in the 5147 A and 5152 A lines, appears at a 
radial velocity of ~ —90 km s^^, which suggests that it 
can be associated with portions of the fast wind seen in 
the most blueshifted P-Cyg absorptions (excluding K I). 
Note also that the velocity of this short-wavelength com- 
ponent is also very close to the velocity of the NAC seen 
in some P-Cyg profiles (see Fig. 15]). As can be seen from 
panel (c) in Fig. [6l the two absorption components of 
the Ti I lines have different shapes. While the blushifted 
component has a very narrow Gaussian profile with a 
FWHM of ~13 km s^^, the redshifted feature is broad, 
with a FWHM of ~40 km s~^, and its profile is asym- 
metric with a steep blue edge and an extended red wing. 

A prominent absorption feature is present at ~7366 A, 
which we have identified as another Ti I line at 7364 A. 
The feature would correspond to the blueshifted com- 
ponent in the discussed above Ti I lines. Note that the 
feature has also a profile similar to the corresponding ab- 
sorptions in the other Ti I lines. A possible redshifted 
component cannot be seen, because it falls in a region 
between the echelle orders. 

We have also identified two weak absorption lines of 
Ca I at 6103 A and 6122 A. They are shown in panel 
(d) of Fig.m These lines are broad (FWHM~60 kms'^) 
and asymmetric. Contrary to the Ti I lines, they have 



an extended blue wing and a sharp red edge. Since these 
Ca I lines arise from a rather high energy level (3.91 eV 
above the ground level), they are expected to be formed 
close to the photosphere of the red supergiant. 

5. MOLECULAR BANDS 

The complexity of the spectrum makes a direct identi- 
fication of molecular bands very difficult. For identifica- 
tion purposes, we have performed a series of simulations 
of molecular spectra. First attempts showed that a model 
capable of reasonably reproducing the observed molecu- 
lar features should include, at least, two components, i.e., 
a stellar photosphere of a cool supergiant and an absorb- 
ing outflowing material with an excitation temperature 
significantly lower than the stellar temperature. 

As a model of the stellar cold photosphere, we have 
used a synthetic spectrum generated from a model atmo- 
sphere from the NextGen grid^ (Hauschil dt ct a l. 1999). 
These models include TiO opacities, what is crucial in 
the case of the analyzed spectrum. We chose a model at- 
mosphere with Teff = 2400K, and with the lowest grav- 
ity available within the grid, i.e. log 5 = 0.0. For these 
parameters there is only a NextGen model with solar 
metallicity, what is probably close to the true metallicity 
of V838 Mon (|KiDDer fc Skodall20?)l . 

Our choice on Tes is a result of initial modeling of 
the observed spectrum, where we used atmosphere mod- 
els for Tcff =2000, 2200, and 2400 K. With the value 
of 2400 K we have obtained most satisfying fits to the 
observed band profiles from higher rotational levels. 

The effective temperature of the cold star can 
also be dete rmined from photometric measurements. 
iMunari et al.l (2007) measured B VRcIc magnitudes o f 
V838 Mon on 2005 Dec. 1, while iHenden et all ()2006D 
measured UBVRJc on 2005 Dec. 25 and JHK' mag- 
nitudes on Dec. 18. Thus the photometry was done 
~2 months after our spectroscopy. Given very slow evo- 
lution of the object at that epoch we can safely assume 
that the magnitudes did not change significantly between 
October and December 2005. Applying the same proce- 
dure to t he ab ove photometric results as that used in 
iTvlendal ()2005n . i.e. fitting a standard supergiant pho- 
tometric spectrum with Eb-v = 0.9 and taking into 
account a contribution from a standard B3V star, we 
have obtained a spectral type of M6-7 and T^s ~ 2480 K. 
Adopting R = 3.1 and a distance of 8 kpc (as in Tvlcnda] 
|2005[) the effective radius and luminosity of V838 Mon 
are ~ 840 i?0 and ^ 2.4 x lO'^L©, respectively. The re- 
sultant fit is shown in Fig. [T] When fitting a blackbody 
to the IJHK magnitudes, a temperature of 2320 K has 
been obtained and, correspondingly, an effective radius of 

905 Rq and a luminosity of 2.2 x 10^ i©. The black- 
body fit (with the contribution from the B3V companion 
added) is shown as a dashed curve in Fig. [T] Thus the 
estimates of Teg obtained from the photometric measure- 
ments are well consistent with the effective temperature 
we used to model the cold photosphere. 

The photosphere model spectrum was then used as an 
underlying spectrum absorbed in the outflowing molec- 
ular gas. In radiative transfer modeling of absorption in 
the outflow (scattering was not included), the gas was 
treated as a plane-parallel homogeneous slab of matter. 

® [ftp : //ftp ■ hs . uni-hamburg ■ de/pub/outgolng/phoenix/NG-giant | 
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Fig. 7. — Model spectra fitted to the photometric 
(UBVRcIcJHK' - full points) measurements of V838 Mon 
in December 2005. Dotted curves - contributions from the 
supergiant and the B3V star. Full curve — final fit (sum of the 
dotted ones). Dashed curve — fit of a blackbody to the IcJHK' 
measurements. See the text for more details. 

In this one-dimensional approach three parameters char- 
acterize the outflowing gas, i.e. column density of an ab- 
sorbing molecule, excitation temperature (the same for 
rotational and vibrational components), and a radial ve- 
locity of the slab (no velocity gradient was implemented). 
Finally, to the resultant spectrum a photospheric contri- 
bution of a B3V star was added (a synthetic spectrum of 
Tefi = 18 000 K and logg = 4.0 was used for this purpose 
- see Sect. [S])- 

It should be noted that the above approach was not 
meant to reproduce the observed spectrum in the whole 
spectral range. Simulations were performed separately 
for different narrow spectral regions and usually for one 
molecule at once. It allowed us to identify and study indi- 
vidual bands dominating in a given spectral range. With 
this modeling procedure rough estimates of physical pa- 
rameters of the outflowing material (excitation tempera- 
ture, column densities, radial velocity) were also possible 
in a number of cases. 

A list of the identified molecular bands is presented in 
Table [3] The first column of the table gives observed 
wavelengths of features. In most cases this is the wave- 
length of a bandhead formed by the R branch, but in 
a few cases this corresponds to a local minimum of the 
absorption feature formed by Q or P branches; otherwise 
an empty space is left. Column (2) contains laboratory 
wavelengths of the features. Columns (3), (4), and (5) 
specify the molecule name, the electronic system with a 
vibrational band identification, and the branch designa- 
tion, respectively. Column (6) gives a heliocentric veloc- 
ity of the feature, if a clear bandhead is observed. Ref- 
erences to appropriate molecular data and comments on 
individual features can be found in the last two columns. 
Most of the identified molecular features are also indi- 
cated on the spectrograms in Figs[lH3l 

The spectrum is dominated by strong TiO absorption 



bands, which are present in all parts of the HIRES spec- 
trum. Prominent are also bands of VO, AlO, ScO, and 
YO. In the following subsections, we describe in some 
detail the identified features for each of the molecules. 

5.1. TiO 

Much effort was devoted to u nderstand elec tronic tran- 
sitions of TiO in labo ratory ([P hillips 1973h. in astro- 
nomical observations (.Ram et al.. .1999.). as well as in 
theoretical computations (|Schwenkelll998f) . In the iden- 
tification procedure, we mainly consulted databases of 
line positions computed from laboratory m easurements 
(|Jorgensenlll99l |PTe3[T99l lSchweny[T99l l. These line 
lists also provide data for bands not measured in labo- 
ratory and in these cases discrepancies between different 
sources are often significant. The lack of accurate spec- 
troscopic data disables a detailed analysis of some TiO 
bands, e.g. the 7' (3,0) band close to 5350 A and higher 
overtones of 7' bands. 

We have identified three electronic systems of TiO in 
the spectrum of V838 Mon, i.e., a (C^A - X^A), 7 
(A^^-X^A), and 7' (B^H-X^A). The molecular bands 
of the 7 and 7' systems have generally a characteristic 
triple structure formed by their main branches. In the 
case of the a transitions, the separation of branches is so 
small that the three components overlap and the triple 
structure is not seen. Additionally, in some bands of the 
7 and 7' systems, relatively strong satellite branches are 
observed, e.g. the R31 branch of the 7' (1,0) band at 
5827 A. Finally, we have also identified forbidden tran- 
sitions of TiO in the spectrum. The latter point is dis- 
cussed in Sect. 15.61 

The observed features of TiO arise both in the stel- 
lar photosphere and in the outflow. Although the latter 
component dominates formation of the observed bands, 
it would be impossible to satisfactorily reproduce a num- 
ber of the observed features without a contribution from 
the photosphere. Unfortunately, we have not been able 
to set valuable constraints on the effective temperature, 
since there are only a few regions in the spectrum, where 
the stellar component dominates. Performed simulations 
show, however, that the selected atmospheric model with 
Tcff = 2400 K works good in reproducing high excitation 
features seen in the spectrum (see below). 

Nonetheless, the most prominent absorptions of TiO 
seen in the spectrum of V838 Mon arise in a low exci- 
tation environment. A comparison between the observa- 
tions and spectra computed from our slab model with dif- 
ferent excitation temperatures suggests a mean value of 
Tex = 500±100K. Adopting this temperature and com- 
paring the observed spectrum to the simulations of the 
7' (2,0), (1,0), (0,0), and 7 (2,0) bands, an estimate of 
the TiO column density of '^16.8 dex cm~^ has been ob- 
tained. Assuming that all Ti is locked in TiO and taking 
a Ti atomic abundance of 10"'' relative to hydrogen, a 
column density of hydrogen of ^23.8 dex cm~^ has been 
derived. 

A characteristic feature of the observed TiO spectrum 
is a strong saturation of bands originating from the 
ground vibrational state. The effect is clearly seen in 
the 7' (0,0) and (1,0) bands. In the 7' (2,0) band the 
saturation is weaker and completely disappears in the 
(3,0) band. Similarly, a strong saturation is seen in the 
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Fig. 8. — Rotational structure of highly excited levels of two 
7 bands of TiO. The observed spectrum (black line) is compared 
with simulations (red line). In order to fit the simulated rotational 
features to those in the observed spectrum, the simulated spectra 
have been shifted to a velocity of V(, = 58kms^^ and smoothed 
with a Gaussian with FWHM = 34 km s~^. An emission feature of 
ScO can easily be seen in the bottom panel (ScO was not included 
in the simulations shown here). 

a (0,0), (1,0), (2,0) bands, and probably also in the (3,0) 
and (4,0) bands. In the red part of the analyzed spec- 
trum, the 7 (0,0) and (1,0) bands are saturated. Due 
to the strong saturation in the 7 and 7' systems the 
triple structure of the bands can be seen only through 
weak emission-like features that appear between the sat- 
urated components, see e.g. the 7 (0,0) and (1,0) bands 
at 6670 A and 7100 A (Fig. [31). 

Radial velocities of some TiO bands can be deter- 
mined by comparing rotational features of the observed 
bands with simulated spectra. Although the rotational 
structure is completely smoothed in the observed band- 
heads, it is recognizable in those parts of the spectrum, 
where only absorptions from higher rotational levels are 
present. Two examples of such rotational features are 
shown in Fig. [51 where a simulated spectrum is over- 
plotted on the observed one. The cold slab had a ex- 
citation temperature of 500 K in the simulations. Note 
however that the contribution from the slab is negligi- 
ble in the 7(3,1) and (4,2) bands. To fit the rotational 
components the modelled spectrum has been shifted to 
= 58±5kms~^. Since these rotational features com- 
ing from high excitation levels are predominantly formed 
in the atmosphere of the central star, the derived velocity 
is expected to be close to the radial velocity of the star. 

Radial velocities can also be determined from the ob- 
served positions of bandheads. In the case of the TiO 
bands present in our spectrum formation of a bandhead 



is often strongly influenced by presence of satellite bands 
and a non-negligible isotopic shift in wavelength. An 
example is presented in panel (a) of Fig. [9l This head, 
belonging to the a (2,0) band, is formed by the R2 branch 
and a contribution from satellite bands is small in this 
case. However, as can be seen in the rotational profile^° 
overplotted on the spectrum in Fig[51 the head is actu- 
ally formed by multiple components. The strongest ones 
belong to the three main TiO isotopomers, i.e., ''^^Ti^^O, 
'i'^Tiif^O, and '^^Ti^^O. For solar Ti isotopic ratios, i.e. 
9.3,10.1,13.4,13.8 for 48/46,48/47,48/49,48/50, respec- 
tively (|Cowievi ll995). all the isotopic components should 
significantly contribute to the shape of the a (2,0) band- 
head. Since they are not distinguishable in the observed 
smooth head, only a position of the most blueshifted 
component belonging to the ^^TiO isotopomer can be 
measured with a reasonable accuracy. In this way a ve- 
locity of Vh = — 67kms~^ has been inferred. 

Another example of a TiO complex bandhead is shown 
in panel (b) of Fig. [9l As seen in the rotational profile, 
the structure of the 7' (1,0) bandhead is even more elab- 
orate than in the previous example. This is due to the 
presence of multiple heads formed by satellite branches. 
The observed head is again very smooth and individ- 
ual components, clearly distinguishable in the overplot- 
ted rotational profile, are not seen in the observed spec- 
trum. An estimate of the outflow velocity is very uncer- 
tain in this case, but by fitting the rotational profile to 
the overall absorption profile of the band, a velocity of 
Vh = —58 ± lOkms"^ has been derived. 

As can be found in Table [31 the most blueshifted 
molecular component is observed in the a (0,0) band 
at 5164.6 A. A velocity of —125 km s^^ has been derived 
assuming that the blue edge of the head is formed by the 
*^TiO isotopomer, instead of the most abundant ^®TiO. 
A similar velocity was found for the a (1,0) band at 
4952. 6A. 

The derived radial velocities of TiO features are in a 
range between 58kms~^ and — 125kms~^. It is reason- 
able to assume that the molecular absorption arises in 
an outflow with a velocity gradient. This can explain 
different velocities of the bands with different excita- 
tion requirements, and can also account for the observed 
smooth structure of the bandheads. 

5.2. AlO 

We based our identification procedure of the AlO 
bands on an analysis of the B ^S"*" - X^S''' system 
performed in iCoxon fc Nax akis* (1985). Absolute 
strengths of the elect ronic transitions ware taken from 
iPartridge et~all (|1983l ). 

The blue-green (B'^'E+ -X'^T,+) band system of AlO is 
evidently present in the analyzed spectrum. The most 
prominent feature is the saturated (0,0) band with a 
main head at 4842 A, see panel (c) in Fig. [9l Three 
bandheads of different branches are close enough to form 
this relatively sharp head. In this case a radial ve- 
locity can be measured accurately and we have found 
Vh = -82±6 km s^^. 

^'^ A rotational profile means here a computed and smoothed 
profile of normalized absorption features of a band and it should not 
be confused with a full simulation from our model involving stellar 
photosphere and an absorbing cool slab as described in Sect. \E\ 
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Fig. 9. — Examples of observed (black line) bandheads of TiO, 
AlO, and VO, together with their simulated rotational profiles 
(red). On each panel a synthetic spectrum of a B3V star is also 
shown (blue). The rotational profiles were shifted in velocity to fit 
the observed bandheads. (a) The a (2,0) bandhead of TiO. The 
three heads clearly distinct in the simulated profile correspond to 
46rpjQ^ 47rpjQ^ iSr^^Q jsotopomcrs. The profile was shifted in 
velocity by -67kms-l. (b) The 7' (1,0) bandhead of TiO. This 
head is formed by multiple strong satellite branches and defining 
a reference laboratory position of the blue edge of the absorption 
feature is problematic. (c) The sharp B-X (1,0) bandhead of 
AlO. In this case the reference laboratory position is well defined. 
The velocity shift is -82kms-l. (d) The C-X (0,0) bandhead 
of VO. As seen in the computed profile, this head has a double 
structure when the band is not very saturated. The velocity shift 
is — 78kms~^. 



Higher vibrational overtones, Av=+1 and Av = +2, 
of the B"X system can also be recognized in the spec- 
trum of V838 Mon. The more intense Av = +1 overtone 
is seen up to the (6,5) band around 4754 A. The radial 
velocity of the bandheads of the first overtone cannot be 
defined accurately due to the poor quality of the spec- 
trum in the blue range. However, the wavelengths of 
bands from higher vibrational levels imply a velocity of 
Vh — 58kms~^. The Av = -|-2 overtone is only hardly 
seen in the spectrum. The most intense (2,0) band is 
contaminated by the He I 4471 A line of the B3V star. 

In the green part of the spectrum, vibrational bands 
belonging to the Av = — 1 sequence are clearly visible. 
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Fig. 10. — Observed spectrum (black line) of the A-X Av=— 1 
subtone of AlO together with its rotational profile computed for 
an excitation temperature of 1500 K (red). The profile was shifted 
in velocity by — 82km~^. The dotted lines mark positions of the 
bandheads if shifted to Vj, = — 82kms~^, while the solid vertical 
lines are drawn for Vh = 58kms~^. The blue line marks the flux 
level of the B3V companion. 



The bands are shown in Fig. [TOl The (0,1) band at 
5079 A is particularly strong. Its bandhead is saturated 
and shifted to Vh = — 82kms~^. One can also distin- 
guish other members of the same overtone originating 
from higher vibrational levels up to v" = 4. According 
to our simulations, these bands, to be observed, require 
a rather high excitation temperature of about 1500 K, 
which suggests that they originate in, or close to, the 
atmosphere of the star. However, sharp bandheads be- 
longing to transitions with v" < 3 are shifted to a velocity 
Vh = — 82kms~^, which is inconsistent with the stellar 
velocity. The velocity rather indicates that the absorbing 
gas is located in outer parts of the outflow. Note that 
this velocity is almost the same as that of the narrow ab- 
sorbing component seen in several atomic lines, e.g. Rb I 
(see Sect. 1121) • 

5.3. VO 

The spectrum of V838 Mon shows strong vibrational 
bands of the C^T,- - X'^T,- and B^H-X^S" systems of 
VO. For the purpose of our simulations a line list for 
these syste nis was compu t ed us ing spectroscopic con- 
stants from iCheung et al.l (119821) Ro tational strength 
factors were taken from iKovacd ( 19690 with corrections 
applied after lWhiting et al.l ( 1975? ). In the case of the B- 
X (0,0), (1,0), and (2,0) bands, a list of energy levels was 
generated by diagonali z ation of t he perturbed Hamilto - 
nian from I Adam et all (|1995l ) and ICheung et all (|1994D . 
The transformation matrices were next used to compute 
rotational strengths of individual lines to properly in- 
clude effects of perturbations. The rotational structure 
of the (2,0) band was not analyzed in laboratory. There- 
fore spectroscopic parameters for this band were extrapo- 
lated from the lower overto nes. Electronic strengths were 
taken from lKarlssonI (|1997f ). Franck-Condon factors were 
computed with the LeRoy's codes" (Le Roy 2004, 2007). 

As can be seen in panel (d) of Fig. [51 the unsaturated 
C-X (0,0) band has two heads. One is formed by the Ri 
and R4 branches, while the other, redshifted by lA, is 
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formed by R2 and R3. Both are present in the observed 
spectrum and the former one appears at a velocity of 
— 77kms~^. A regular ladder of rotational lines of the 
(0,0) band formed by main P branches can be recognized 
starting at 5753 A (see Fig. [5]). Fortunately, precise lab- 
oratory wavelengths are available for these lines and a 
detailed modeling of the band was possible. The simu- 
lations suggest that this part of the spectrum is formed 
at a low excitation temperature of '^SOOK. The derived 
velocity of the absorbing gas in the rotational transitions 
of the (0,0) band is V/j =43kms~^. From the heads of 
higher transitions of the C-X system, namely the (1,0) 
and (2,0) bands, velocities of —77km s~^ and — 45kms~^ 
were derived, respectively. 

Additionally, within the B~X system, the bands (0,0), 
(1,0), and (2,0) were identified. The (0,0) and (1,0) 
bands, together with their simulated spectra, are pre- 
sented in Fig. [11] Most of the features of these bands are 
only slightly shifted from their laboratory wavelengths 
and the measured velocity is V;i = — 2kms~^. A rather 
uncertain value of ^^BSOK was deduced from the sim- 
ulations as an excitation temperature of the absorbing 
gas (mainly from the intensities of absorption features 
at 7910 A and 7920 A). Furthermore, in order to repro- 
duce some features seen in the spectra of the bands with 
higher vibrational numbers, e.g. the absorption feature 
at 7474A belonging to the (2,1) band, it was necessary 
to add an additional layer of absorbing VO gas with a 
rather high excitation temperature, namely Tcx=1500K. 
This hotter component, treated in the simulations as a 
slab between the stellar photosphere and the cold slab, 
appears at a velocity of ^^SSknis^^. 

Our simulations of the C-X (0,0), B-X (0,0) and (1,0) 
bands give satisfactory fits to the observations when a 
column density of the VO absorber is of ~16.4 dex cm~^. 
If we assume that all V is locked in VO molecules and tak- 
ing the solar abundance of V relative to H of 10~^, this 
corresponds to a column density of H of ^24.4 dex cm~^. 
The result is consistent with a similar estimate performed 
from the TiO bands in Sect. 15. ll 

5.4. ScO 

The only electronic system of ScO present in the spec- 
trum is the A^n — X^I]+ system. Absorption features of 
the (0,0), (1,1), (2,1) bands are clearly seen in the green 
part of the spectrum. Moreover, one band of ScO, i.e. 
the (0,1) band, was found in emission. For the purpose of 
our simulations, line positions of the A-X (0,0) band were 
computed from the Hamiltonian diagonali zation parame - 
terized bv lStringat et all (|1972i ) and iAdams et all (|1968D . 
Higher vibrational levels were approximated from data 
in NI ST^^. For the positi ons of transitions of the (1,0) 
band, IRice fc Field (|T986l) was consulted. Strengths of 
rotational lin es were found in the approach proposed in 
iHerbigj ()1974[ ). It should be noted, that although the line 
positions of the (0,0) band are exact, their strengths are 
only approximate. 

Figure [12] shows a region of the observed spectrum cov- 
ering two subbands of the (0,0) band. Simulations of this 
spectral region suggest the presence of at least two differ- 
ent velocity components. The bandhead formed by the 

1^ |http : //physics ■ nlst . gov/PhysRef Dat/ASD/ 1 
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Fig. 11. — Comparison between the observed (black line) and 
simulated spectra of the B-X bands of VO. Top: the (0,0) band 
was simulated with the excitation temperatures of 350 K (red 
line) and 500 K (green line). The simulated spectra were shifted 
to a velocity of -2kms~'^ and smoothed with a Gaussian with 
FWHM = 65 km s"'^. The vertical lines mark positions of the band- 
heads for a radial velocity of -80 km s~^, which fits to the observed 
bandheads of ^Q42 and R2, and to the blue edge of the Rb I pro- 
file. Bottom: the simulation of the (1,0) band obtained with two 
layers of the absorbing molecular gas with different temperatures, 
i.e. 500 K and 1500 K. The colder layer was shifted in velocity to 
— 2kms~^, while the warmer one, treated in our radiative trans- 
fer modeling as a layer between the stellar atmosphere and the 
cold slab, was shifted to a velocity of 58kms~^. The dotted lines 
mark positions of the bandheads for a radial velocity of —60 km s ~ ^ , 
which fits to the ^R42 and ^Qsi bandheads. The absorption fea- 
ture at 7474 A belongs to the (2,1) band. 



'n3/2 ^Q2G + ^i?2G (at 6036 A) and ^Hi/a ^Qia + '^RiG 
(at 6079 A) branches are shifted to — 55kms~^ (as mea- 
sured for the lowest part of the profile). The second 
velocity component at 45kms~^ was found from the po- 
sition of the head formed by the ^RiG branch at 
6064. 17 A. This latter head is formed at a rather high 
rotational number, i.e., N=58, contrary to the heads ob- 
served at — 55kms~^, which are formed at N=28 and 
N=17. The relative velocity shift can be understood as 
a result of decreasing population of high rotational levels 
with increasing velocity in the observed outflow. A lack 
of accurate laboratory wavelengths disables an analysis 
of other absorption bands of ScO. 

Remarkable is the presence of the A-X (0,1) 
^1^3/2 ^Q2G and ^ni/2^QiG subbands of ScO in emis- 
sion at 6408 A and 6458 A, respectively (see bottom panel 
of Fig. [5] for the latter subband). The observed veloc- 
ity of this features agrees well with the velocity of emis- 
sion components seen in the atomic P-Cyg profiles (see 
Sect. 14. 2p . The origin of this molecular emission can be 
explained by fluorescent pumping of the upper electronic 
level by absorption in the (0,0) 6090 A and 6132 A sub- 
bands, clearly seen in the spectrum of V838 Mon. 
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Fig. 12. — Comparison between the observed and simulated spec- 
trum of the A— X (1,0) band of ScO. The two upper plots are rota- 
tional profiles of the band, computed for excitation temperatures 
of 750 K and 500 K (blue and red, respectively). They are shifted 
to a velocity of 58kms~^. The observed spectrum is drawn with 
a black line. Overlaid on it in red is a simulated spectrum with 
a molecular slab absorbing the light of the cool star as described 
in Sect. |5] The absorbing layer has an excitation temperature of 
500 K and is shifted in velocity by -55 km s~^ . The simulated spec- 
trum is smoothed with a Gaussian with FWHM = 32 km s~^. The 
flux level of the B3V star is indicated with the blue line. 



5.5. YO 

All the features of YO identified in the spectrum 
of V838 Mon belong to the A'^n-X'^j:+ system. A 
list of the transitions used in our identification pro- 
ce dure was prepared o n the basis of data provided 
in [Bernard fc Gravinal (|1983f ). There is another well 
known system of YO , i.e. B^II - X'^E+ around 4819 A 
(jBernard et al.l[T979[ l. but it falls into a spectral region of 
saturated TiO bands and was not found in our spectrum. 

We performed a detailed simulation of the A^n3/2 — 
X^I]+ transition. The results of this simulations are 
shown in Fig. [TS] The upper two plots in Fig. [13] are 
rotational profiles of the transition obtained for two ex- 
citation temperatures, 750 K and 500 K. It can be seen 
that a sharp edge of the (0,0) bandhead at 5973 A is 
formed with the higher temperature. In the observed 
spectrum this head is very smooth and resembles more 
the feature seen in the profile for 500 K, and this is a 
value we took for further simulations. It should be noted, 
however, that the observed profile of the band can be in- 
fluenced by the physical structure of the wind, which is 
ignored in our simulations; in particular, a velocity gra- 
dient in the molecular outflow can smooth the shape of 
the bandheads, and even if the temperature is high, no 
sharp head would be observed. Nonetheless, to perform 
full simulations of the A^n3/2 — X^I]+ transition of YO 
it was necessary to include the 7 (3,0) band of TiO as 
an additional important absorber in this spectral range. 
The green line in Fig. [13] shows the TiO contribution 
to the observed spectrum. The final simulation which 




5970 5980 

Wavelength (A) 



Fig. 13. — Same as in Fig. 1121 but for the fine components of 
the electronic transitions of the A^n3/2 — X'^S^' feature of YO. 
The final simulation (bottom) includes an extra absorption of TiO, 
which is shown with a green line. The simulation was performed 
with a slab of 500 K and the stellar photosphere. All the computed 
spectra were shifted to a velocity of 58kms~^ and smoothed with 
a Gaussian with FWHM = 37 km s-^. 



includes absorption of YO and TiO with an excitation 
temperature of 500 K is shown at the bottom of Fig. [13] 
From the overall fit of the (0,0) band to the full simula- 
tion a velocity of Vh =-8kms~^ has been derived. 

Note however that a sharp (1,1) bandhead formed by 
the ^Q2i branch at 5988 A appears at higher excita- 
tion temperatures, at least 750 K, as can be seen from 
Fig. [T2] Such a sharp feature seems to be present in 
the observed spectrum, suggesting that this absorption 
feature is formed in gas with a somewhat higher temper- 
ature than the one found for the (0,0) bandhead. This 
is confirmed by the measured velocity of the (1,1) band- 
head, Vfi — 58kms~^, which is the same as that of other 
high-excitation features found in the spectrum. 

5.6. TiO - forbidden transitions 

Very interesting from the spectroscopical point of view 
is the presence of forbidden transitions of TiO, i.e. the 
ci$-x3A (0,0) band around 4737A and the b^H-X^A 
(0,0) band at 6832 A. The stronger b-X band is shown 
in detail in Fig. [TJ] The intensities of the bands were 
predicted in ab initio calculations of ro-vibrational en- 
ergy level s of th e 13 lowest electronic states of TiO by 
iSchwenkd (|1998D . As shown in Fig. [H our simulation 
can reproduce the forbidden bands quite well. The best 
fit indicates a velocity of Vh — — 80kms^^ and an exci- 
tation temperature of about 500 K. However, when the 
column density found for other TiO bands (see Sect. 15. Ij ) 
is assumed, oscillator strengths of the b^II— X^A band 
must be enhanced by a factor of 70±35 in order to re- 
produce the observations. To our knowledge this is the 
first, in astrophysics as well as in laboratory, observation 
of TiO forbidden bands. 
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Fig. 14. — Observed and simulated spectrum of the forbidden 
band b^H-X^A (0,0) of TiO. The simulation (red line) was ob- 
tained with an excitation temperature of 500 K. It was shifted to 
a velocity Vh=-82kms~^ and smoothed with a Gaussian with 
FWHM = 32kms-i. 



6. PHOTOSPHERIC SPECTRUM OF THE B3V 
COMPANION 

A number of broad absorption lines of the Balmer se- 
ries, from H/3 down to at least Hl2, as well as a few He I 
lines, can clearly be seen in the blue part of our spectrum 
(Fig.[T|). These features are attributed to the B3V com- 
panion. We compared the observed spectrum to a grid 
of synthetic sp ectra BSTAR2006, gene rated with the Syn- 
spec program (|Lanz &: Hubenvl 120071 ) . Solar metallicity 
models, with effective temperatures between 15 000 and 
30 000 K and gravities in the range logg — 2.0 — 4.75 
were considered. A turbulent velocity of 2 km s~^ was 
adopted. Good fits to the observed spectrum were ob- 
tained for Toff = 18 000 ± 2 000 K and log 5 = 4 ± 0.75. 
The low signal-to-noise ratio in the part dominated by 
the B3V star and uncertainties in the flux calibration did 
not allow us to obtain more precise estimates of the star 
parameters. In any case, the above values agree well with 
those expected for a B3V star. In the following discus- 
sion we use the synthetic spectrum with T^s = 18 000 K 
and \ogg = 4.0. This spectrum, broadened with a ro- 
tational velocity of 250 km s~^ (see below) and shifted 
to a radial velocity found for the hot star (see below), is 
overplotted on the observed spectrogram of V838 Mon 
in Figs. [DEI 

A comparison of the observed spectrum to the syn- 
thetic one shows that the absorption lines of the B3V 
companion are rotationally broadened. In order to find 
the rotational velocity, we convolved the synthetic pro- 
files with a Gaussian profile broadened with a grid of ro- 
tational velocities (we made use of the roim5 program^^). 
Note that this procedure does not take into account limb 
darkening nor gravitational darkening, which can be im- 
portant for profile formation in a rapidly rotating star. 
Lines usually used to investigate rotational broadening 
of B-type stars, i.e. He I 4471 A and Mg II 4481 A, are 
weak, noisy, and affected by emission lines and/or molec- 
ular bands in our spectrum. However, we have found that 
in the green part of the spectrum the He I 5876 A line 

htt p : //nova ■ astro ■ umd ■ edu/Synspec43/ | 
sjrnspec-f rames-rotin . html 



is relatively strong and has a "clean" profile. This green 
spectral region is dominated by the cool photosphere but, 
fortunately, in the spectral vicinity of the A5876 line the 
flux from this component is practically reduced to zero 
due to a saturated TiO band (see Fig. [2|) and the line 
from the B3V component is well seen. Thus, when deter- 
mining the rotational velocity of the B3V component we 
relied mainly (but not only) on the profile of He I 5876 A. 
Using the Toff =18 000 K, log 3=4.0 template we found a 
projected rotational velocity V sn\i — 250 ± 50 km s~^. 
This is a rather large value^** and indicates that the ro- 
tation axis is close to the plane of the sky. 

We also attempted to measure the radial velocity of 
the hot star. We used the Fourier cross-correlation 
method (jTonrv fc David Il979( ) as it is implemented in 
the IRAF's task rvcor. As a template, we used our syn- 
thetic spectrum broadened to the rotation velocity of 
250 km s~^. The cross-correlation was performed for 
the region between 3815 A and 4180 A, with all emis- 
sions and ISM absorptions in this range removed man- 
ually. (We performed tests with ranges including H/3 
and H7, but they did not give reliable results. These 
spectral regions are too strongly contaminated by atomic 
emissions and molecular absorptions. Also, including 
the very noisy shortward part of the spectrum spoils 
the cross-correlation function and gives unrealistic re- 
sults.) As a result we obtained a heliocentric velocity 
of 37.0 ± 16.8 km s~'^ [with a Tonry & Davis' (1979) r 
value of 18.0], or correspondingly 19.6 ± 16.8 km s^^ in 
to the local standard of rest. The error value was taken 
directly from the output of the rvcor procedure and, in 
our opinion, it is slightly overestimated. 

7. INTERSTELLAR FEATURES 

Several interstellar absorption features can be identi- 
fied in our spectrum. The strongest ones are those of 
Ca II H, K and Na I Di, D2. They are shown in Fig. [I5l 
Interstellar components can also be found in the deep 
P-Cyg absorption of K I at 7699 A (see upper panel in 
Fig. EI). The Na I and K I interstellar absorptions in the 
spectrum of V838 Mon hav e already been studied in a 
numb er of papers (e.g. fKipper et al.l 120041 : IMunari et all 
[2OO5I) . As can be seen from Fig. 1151 the profiles of the 
Ca II lines, not described in the literature so far, have 
a similar structure to those of Na I. All the four profiles 
displayed in Fig. [15] show three compone nts at 32, 42, 
and 65 km s~^, similarly as observed in IMunari et al.l 
([2005 ). We note, however, that unhke previous studies, 
where the interstellar lines were observed in the spectrum 
of V838 Mon itself, the features presented in Fig. [TH] are 
observed in the spectrum of the B3V companion. As can 
be seen from Figs. [1] and [2] the local continuum near the 
discussed lines comes from the B3V companion. At the 
wavelengths of the Ca II lines V838 Mon is too cool to 
emit any significant radiation, while at the wavelengths 
of the Na I lines an absorption band of TiO is so strong 
that it eats up practically all the flux from the cool su- 
pergiant. The fact that the interstellar line profiles in 
the spectrum of the B3V companion are very similar to 

For instance, in Abt et al. (2002) one finds that the mean 
projected rotational velocity for 106 dwarfs with spectral types 
B3-B5 is 108±8 km s'^, and only 2.4% of the stars have V sini > 
300 km s-i. 
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Fig. 15. — Interstellar doublets of Ca II and Na I. The curves, 
except Na I Di, are shifted in the vertical scale to better show 
the profile structures. The Ca II H profile is noisy as the line is 
observed near the core of the He absorption, where the stellar flux 
is very low (see Fig.[T]l. Three main absorption components at 32, 
42, and 65 km s~^ are marked with dashed vertical lines. 

those in the spectrum of V838 Mon itself is consistent 
with the idea that both stars are related, at least that 
they are at a similar distance. 

We have identified six diffuse interstellar bands (DIBs) 
in our spectrum, i.e., at AA4428, 5780, 5797, 628 4, 6614 
and 72 24 [the central wavelengths are taken from lHerbid 
(fl995l) ]. They aU are redshifted to a heliocentric velocity 
of ~60 km s~^. It is well known that equivalent widths 
of the features at 5780 and 5797 A correlate very well 
with interste l lar ex tinction. Using conversion formulae 
from iHerbigl (fT975l ) we derived Eb-v of 0.95 and 0.90 
for the DIBs at 5780 A and 5 797 A , respectively. Tak- 
ing the relation from Herbid (|1993f) . which is based on 
a larger number of data points, gives Eb-v — 1-00 for 
both features. Given considerable uncertainty of these 
estimates, due to uncertainties in the equivalent width 
measurements and the intrinsic scatter in the relations 
used, our result is consistent with E b-v — 0-9 ob tained 
in other studies [as summarized in .Tvlendal (|2005D ] . 

8. SUMMARY 

We have obtained a spectrum of V838 Mon on 2005 Oc- 
tober 13 using the Keck/HIRES instrument. The spec- 
trum of a resolving power of ~ 34 000 covered a wave- 
length range of 3720 - 7960 A. The spectrum is very com- 
plex and shows numerous spectroscopic features which 
can be identified as arising from different components 
or environments. V838 Mon itself, i.e. the object which 
erupted in 2002, is seen mainly in the green and red parts 
of the spectrum as a cool continuum with strong, broad 
and complex molecular absorptions. Numerous atomic 
lines showing P-Cyg profiles or broad absorption features 
are also present in this component. The blue part of the 
spectrum is dominated by a photospheric spectrum of 



the B3V companion, showing broad absorptions lines of 
H I and He I. In the blue and green parts of the spec- 
trum, numerous emission lines can be identified. These 
mostly arise from forbidden transitions in singly ionized 
iron. Finally, several absorption features of interstellar 
origin have been identified. These include atomic reso- 
nance lines as well as DIBs. 

The numerous [Fe II] emission lines observed in our 
spectrum evidently arise from rarefied and (partly) ion- 
ized matter. All the lines show the same profile, which 
can be well fitted with a Lorentzian profile. The lines are 
centred at a heliocentric radial velocity of ~13 kms~^ 
and have a typical FWHM of -80 kms'^ 

In the cores of the photospheric absorption lines of the 
B3V companion we see weak emission features. The dou- 
ble feature in the H/3 line can be interpreted as a single 
emission with a narrow absorption superimposed on it. 
The derived parameters of the emission component are 
close to those of the [Fe II] lines. This suggests that the 
Balmer emission features arise from the same environ- 
ment as the [Fe II] emission lines. 

The spectrum displays a number of strong atomic, 
mostly resonance, lines with P-Cyg profiles. Evidently 
V838 Mon, being almost 4 years after its outburst, con- 
tinues losing mass. The emission component typically 
peaks at a heliocentric velocity of '^80 kms~^, while the 
absorption component usually extends from —50 kms""'^ 
to — —100 kms^^. The absorption component often 
shows a complex profile. In particular, in a few lines a 
strong and narrow absorption component is clearly seen 
at a velocity of ~ —80 kms~^. A few atomic lines were 
observed to show pure absorption profiles. 

The spectrum from the cool component is dominated 
by numerous, often very strong, molecular bands in ab- 
sorption. Some bands are so deep that the flux from the 
cool component practically goes to zero. In these cases 
we have a sort of spectral windows, through which we 
can see the B3V companion, even in the red part of the 
spectrum. All the identifled bands are from oxides. They 
include TiO, VO, AlO, ScO, and YO. 

Our analysis shows that the bands are formed in differ- 
ent regions. Some of them, especially those from highly 
excited states, are formed, at least partially, in the pho- 
tospheric regions. The temperature in these layers is as 
high as —2400 K. These bands show a heliocentric ra- 
dial velocity of —58 kms^^, which we identify as a radial 
velocity of V838 Mon itself. 

Most of the bands are formed in regions significantly 
cooler than the photospheric ones. The vibrational tem- 
peratures derived from some bands go down to —500 K. 
These are usually very strong bands arising from ground 
states. The radial velocity of these bands is usually much 
lower than that of the photospheric regions and goes 
down to — —125 kms""'^. Overall, our analysis of the 
molecular bands shows that they are mostly formed in 
outflowing matter with a significant velocity gradient. 

We have detected bands arising from forbidden tran- 
sitions in TiO. As far as we know, these are first ob- 
servations of TiO forbidden transitions in astrophysical 
objects. 

The observed spectrum of the B-type companion, seen 
primarily in the blue part of our spectrum, can be fit- 
ted with a model atmosphere spectrum with solar abun- 
dances, effective temperature of —18 000 K, and a grav- 
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ity, log 5 ~ 4.0. These values agree well with those ex- 
pected for a B3V star. The absorption lines are sig- 
nificantly broadened, which we interpret as a rotational 
broadening with a velocity, Vs'mi ~ 250 kms~^. We 
have derived a heliocentric radial velocity of the B3V 
companion of 37 ± 17 kms^^. 

Several interstellar absorption features are seen in our 
spectrum. They include lines of Ca II, Na I, and K I, 
as well as DIBs. The lines of Ca II and Na I are ob- 
served in the spectrum of the B3V companion. Their 
profiles are very similar to those observed in the spec- 
trum of V838 Mon during the 2002 outburst. This result 
is consistent with the idea that V838 Mon and the B3V 
companion are related objects. From equivalent widths 
of diffuse interstellar bands we have estimated an inter- 
stellar reddening oi Eb-v — 0.9 — 1.0, which agrees well 
with interstellar extinction values obtained by other au- 
thors using other methods. 

A detailed analysis and interpretation of the data pre- 
sented in the present paper, as well as of data from other 
sources, will be done in a forthcoming paper. 
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Keck/HIRES spectrum of V838 Mon 



TABLE 1 

Identified emission lines 



Kib (^) Multiplet Flux Notes 



3883.820 [Fc II 


? 


24F 




a 


3993.059 Ni II 


7 


4F 






4101.734 


m 




1 






4177.196 


Fell 




21F 


3.355C-15 


b 


4201.172 [Nill 


7 


3F 






4243.969 


Fell] 


21F 


5.028C-14 


c,d 


4244.813 


Fell] 


21F 


1.448C-14 


c,d 


4276.829 


Fell] 


21F 


5.469C-14 




4287.394 


Fell] 


7F 


1.423e-13 




4305.890 


Fell] 


21F 


1.292e-14 




4319.619 


Fell] 


21F 


2.852«^14 




4326.237 [Ni II 


7 


3F 






4340.464 


H7 




1 






4346.852 


Fell 




21F 




c 


4352.778 


Fell 




21F 


2.448^14 




4358.360 


Fell 




21F 


2.8426-14 


d 


4359.333 


Fell 




7F 


1.022e-13 


d 


4o7z.427 


Fell 




21F 


2.9406-14 




4ooz. r 4z 


Fell 




6F 






441o.7ol 


Fell 




7F 


1.1096-13 


c,d 


441o.2oD 


Fell 




6F 


5.5796-14 


c,d 


/f /( on A AT 

44oz.447 


Fell 




6F 






44oz.uyo 


Fell 




7F 


3.2216-14 




44o / .y4o 


Fell 




6F 


2.6276-14 




/I /I T /I f\r\ A 

4474.904 


Fell 




7F 


1.817e-14 




4450.749 


Fell 




6F 


8.752C-15 




4492. oo4 


Fell 




6F 


4.873e-15 


b 


A crnn i^no 

4ouy.ouz 


Fell 




6F 


7.680e-15 


e 


4514. 9UU 


Fell 




6F 


6.475C-15 


e 


4639.667 


Fell 




4F 


1.049C-14 


f 


4dd4.44u 


Fell 




4F 


3.666C-15 


e 


4728.068 


Fell 




4F 


2.069C-14 




4774.718 


Fell 




20F 


1.120e-14 




,10 1/1 cr o /I 

4ol4.oo4 


Fell 




20F 


4.4956-14 




4obl.ozo 


H/3 




1 






4o74.4oo 


Fell] 


20F 


1.4596-14 




inter- chip gap 












1 o.ooo 


Fell] 


20F 


1.8366-14 


g 


O U (J (J . J- 


Fell] 


20F 


2.2256-14 




5018.440 


Fell 




42 


2.5706-14 


c,d 


5020.233 


Fell] 


20F 


1.1636-14 


c,d 


5041.024 


Si II? 


5 








Fell] 


20F 






5055.984 


Si II? 


5 




c 


5056.317 


Si II? 


5 




c 


5158.001 


Fell] 


18F 


3.7230-14 


c,d 


5158.777 


Fell] 


19F 


2.433e-14 


c,d 


5163.951 


Fell] 


35F 






5169.033 


Fell 




42 


2.411C-14 




5181.948 


Fell] 


18F 


6.598C-15 




5184.788 


Fell] 


19F 


1.915C-15 


b 


5197.902 


[N I] 




IF 




h 


5220.059 


Fell] 


19F 


1.203e-14 


g 


5261.621 


Fell] 


19F 


2.979e-14 


f 


5268.874 


Fell] 


18F 






5273.346 


Fell] 


18F 


5.1606-14 


f 


5333.646 


Fell] 


19F 


2.5066-14 


f 


5347.653 


Fell] 


18F 




c 


5376.452 


Fell] 


19F 


2.388^14 




5477.241 


Fell] 


34F 


8.6806-15 


g 


5527.609 


[Fe II] 


34F 


1.6776-14 


g 


5577.339 


[O I] 


3F 






5746.966 


[Fe II] 


34F 


1.3296-14 




6300.304 


[OI] 


IF 






6363.776 


[OI] 


IF 






inter-chip gap 












7155.160 


[Fe II] 


14F 


3.861c 15 





Note. — The integrated line fluxes were measured 
on dcrcddened spectra and they are given in units of 
erg crri^^ s^^. 

multiplet 24 absent in NIST 
^ weak line; flux uncertain 

blend 

flux measurement obtained with a deblending procedure 



Kaminski et al. 



TABLE 2 
P-Cyg and absorption lines 



Kl (A) 


Ion 


Multiplct 


Profile type 


Notes 


ovoo.bbo 


L-a 11 Iv 


1 


absorption 


a 




Ca II H 


I 






4254.332 


Cr I 


1 


absorption 


b 


4274.796 


Cr I 


1 


absorption 


b 


4289.716 


Cr I 


1 


absorption 


b 


4545.945 


Cr I 


10 


absorption 




4554.033 


Ba II? 


1 


absorption 


c 


4571.096 


Mgl 


1 


P-Cygni 


d 


4580.043 


Cr I 


10 


absorption 




4613.357 


Cr I 


21 


absorption 




4616.120 


Cr I 


21 


absorption 




4646.148 


Cr I 


21 


absorption 


8 


4651.282 


Cr I 


21 


absorption 


8 



inter-chip gap 



5060.078 


Fe I 


1 


absorption 




01 lU.^iO 


Fo I 
re i 


1 
± 


P-Cygni 


f 
I 


5127.681 


Fe I? 


1 


absorption 




5147.479 


Ti I 


4 


absorption 


g 


5152.185 


Ti I 


4 


absorption 


g 


5204.506 


Cr I 


7 


absorption 


e 


5206.038 


Cr I 


7 


absorption 


e 


5208.419 


Cr I 


7 


P-Cygni? 


e 


5247.049 


Fe I? 


1 


absorption 




5264.157 


Cr I 


18 


absorption 


e 


5265.723 


Cr I 


18 


absorption 


e 


5296.691 


Cr I 


18 


absorption 


e 


5298.277 


Cr I 


18 


P-Cygni 


e 


5300.744 


Cr I 


18 


absorption 




5345.801 


Cr I 


18 


absorption 


e 


5348.312 


Cr I 


18 


absorption 


e 


5394.676 


Mn I 


1 


P-Cygni 




5409.772 


Cr I 


18 


P-Cygni 




5426.237 


Ti I 


3 


absorption 


g? 


5432.546 


Mn I 


1 


P-Cygni 




5535.481 


Ba I 


2 


P-Cygni 




5698.520 


V I 


35 


absorption 




5703.580 


V I 


35 


absorption 




5706.980 


V I 


35 


absorption 




5727.050 


V I 


35 


absorption 


e 


5727.650 


V I 


35 


absorption 


e 


5889.951 


Na I D2 


1 


absorption 




5895.924 


Na I Di 


1 


absorption 




6102.723 


Ca I 


3 


absorption 




6122.217 


Ca I 


3 


absorption 




inter- chip gap 










6551.676 


Fe I 


13 


absorption 




6556.062 


Ti I? 


102 


absorption 




6572.779 


Ca I 


1 


P-Cygni 




6707.761 


Li I 


1 


absorption 


e, 


6707.912 


Li I 


1 


absorption 


e, 


7364.097 


Ti I? 


97 


absorption 




7698.974 


K I 


1 


P-Cygni 




7800.268 


Rb I 


1 


P-Cygni 




7947.603 


Rb I 


1 


P-Cygni? 





Note. — Photospheric signatures of the B3V companion 
are not included. 

ISM feature 
^ weak line; flux uncertain 

Cs I 4555.8 A resonance lino? 
absorption component very weak 
" blend 

f emission component weak 
8 two separated components 
^ P-Cygni profile? 



Keck/HIRES spectrum of V838 Mon 

TABLE 3 

Identified molecular bands 





(A) 


Molecule 


Electronic 
transition 




Feature 


Vel. 
(kms-i) 


Ref. 


Notes 


4420.5 


4421.521 


TiO 


a 


(4,0) 


R2 




1 




4469.4 


4470.54 


AlO 




(2,0) 


R 




2 




4494.8 


4494.04 


AlO 


b2e+ - x2s+ 


(3,1) 


R 




2 




4517.5 


4516.40 


AlO 


b2s+ - x2s+ 


(4,2) 


R 




2 




4582.5 


4584.056 


TiO 


a 


(3,0) 


R2 




1 




4624.8 


4626.08 


TiO 


a 


(4,1) 


R2 




1 




4647.7 


4648.23 


AlO 


B - X ZjT^ 


(1,0) 


R 




2 




4671.7 


4672.02 


AlO 


B^S+ - X^S+ 


(2,1) 


R 




2 




4694.5 


4694.62 


AlO 


B^S+ - X2S+ 


(3,2) 


R 




2 




4715.4 


4715.54 


AlO 


B2e+ - X2S+ 


(4,3) 


R 




2 




4736.0 


4735.82 


AlO 


B2e+ - X2S+ 


(5,4) 


R 




2 




4738.0 


4737.460 


TiO 


ci$ - X^A 


(0,0) 


Rl3 




2 




4754.4 


4754.29 


AlO 


B2s+ - X2S+ 


(6,5) 


R 


+58 ± 10 


2 




4757.9 


4759.005 


46TJ160 


a 


(2,0) 


R2 


-67 ± 6 


3 


a 


4759.1 


4759.973 


47TJ160 


a 


(2,0) 


R2 




3 






4760.902 


«Tii«0 


a 


(2,0) 


R2 




3 




4803.0 


4804.333 


TiO 


a 


(3,1) 


R2 




1 


b 


4840.95 


4842.27 


AlO 


B2e+ - X2S+ 


(0,0) 


R2 


-82 ± 6 


2 




4952.6 


4954.562 


TiO 


a 


(1:0) 


R2 


-119 ± 6 


1 




4997.9 


4999.129 


TiO 


a 


(2,1) 


R2 




1 




5044.9 


5044.501 


TiO 


a 


(3,2) 


R2 




4 




5078.0 


5079.36 


AlO 


B2e+ - X2S+ 


(0,1) 


R 


-82 ± 6 


2 




5100.7 


5102.13 


AlO 


B2e+ -X2S+ 


(1,2) 


R 


-82 ± 6 


2 




5122.0 


5123.33 


AlO 


B2e+ -X2S+ 


(2,3) 


R 


-82 ± 6 


5 




5143.7 


5142.89 


AlO 


B2e+ -X2S+ 


(3,4) 


R 


+48 ± 12 


2 




5161.8 


5160.98 


AlO 


B2e+ - X2S+ 


(4,5) 


R 


+48 ± 12 


2 




5164.6 


5166.753 


TiO 


a 


(0,0) 


R2 


-125 ± 12 


1 


c 


5227.4 


5228.2 


VO 


CS- - X^E- 


(2,0) 


R4 


-45 ± 12 


2 




5276 


5275.8 


VO 


C-*S- - X^E- 


(3,1) 






2 




5308.4 


5307.337 


TiO 


a 


(3,3) 


R2 




4 




5312 


5310.928 


TiO 


a 


(3,3) 


R3 




4 




5319.0 


5319.746 


TiO 


7' 


(3,0) 


Ri 




4 


d 


5320.4 


5321.111 


TiO 


7' 


(3,0) 


Qi 




4 


d 


5338.6 


5341.240 


TiO 


7' 


(3,0) 


R2 




4 


d 


5340.4 


5343.223 


TiO 


7' 


(3,0) 


Q2 




4 


d 


5356.5 


5355.975 


TiO 


a 


(4,4) 


R2 




4 




5361.5 


5363.528 


TiO 


7' 


(3,0) 


R3 




4 


d 


5365.0 


5367.167 


TiO 


7' 


(3,0) 


Qs 




4 


d 




5376.83 


AlO 


B2s+ - X2S+ 


(2,4) 


R 




2 






5411.960 


TiO 


7' 


(4,1) 


R3 




4 






5415.639 


TiO 


7' 


(4,1) 


Qs 




4 




5447.2 


5448.233 


TiO 


a 


(0,1) 


R2 


-57 ± 18 


1 




5467.9 


5469.3 


VO 


C^E- - X^S- 


(1,0) 


R4 


-77 ± 12 


2 




5497.7 


5496.742 


TiO 


a 


(1,2) 


R2 


+55 ± 12 


1 




5500.4 


5499.741 


TiO 


a 


(1,2) 


R3 


+36 ± 12 


1 






5517.3 


VO 


C-*E- - X^S- 


(2,1) 


R4 




2 






5546.041 


TiO 


a 


(2,3) 


R2 




1 




5562.1 


5562.130 


TiO 


7' 


(2,0) 


■5R21 




5 




5568.9 


5569.089 


TiO 


7' 


(2,0) 


Ri 




5 




5591.4 


5592.609 


TiO 


7' 


(2,0) 


R2 




5 




5616.4 


5616.792 


TiO 


7' 


(2,0) 


Rs 




5 




5640.3 


5640.849 


TiO 


7' 


(3,1) 


R2 




5 




5662.0 


5662.223 


TiO 


/3 


(2,2) 


R 




5 




5662.8, 5665.0 


5665.19 


TiO 


7' 


(3,1) 


Ra 




5 




5669.5 


5669.966 


TiO 


7' 


(3,1) 


Qa 




5 




5688.5 


5689.640 


TiO 


7' 


(4,2) 


R2 




5 




5691.3 


5692.270 


TiO 


7' 


(4,2) 


Q2 




5 




5713.8 


5714.805 


TiO 


7' 


(4,2) 


Ra 




5 




5718.4 


5719.223 


TiO 


7' 


(4,2) 


Qa 




5 




5735.2 


5736.703 


VO 


C-*E- - X^E- 


(0,0) 


R4 


-78 ± 6 


2 




0/09.4, 5/00.4 




TiO 


a 


(0,2) 


-K.2 




3 




5773.2 


5772.74 


ScO 


A^Hg/a - X2s+ 


(2,1) 


«Q21 




2 


e 




5811.60 


ScO 


A^Hi/j - X2S+ 


(2,1) 


Qi 




2 


e 


5811.0 


5809.872 


TiO 


a 


(1-3) 


R2 




3 




5827.7 


5827.024 


TiO 


7' 


(1:0) 


Rai 




3 




5837.3 


5839.005 


TiO 


7' 


(1:0) 


R21 


-58 ± 10 


3 






5847.593 


TiO 


7' 


(1:0) 


Ri 




3 




5948.2 


5947.680 


TiO 


7' 


(2,1) 


Ra 




3 




5956.2 


5954.419 


TiO 


7' 


(2,1) 


Qa 




3 




5972.5 


5972.17 


YO 


A2n3/2 - x2e+ 


(0,0) 


«Q21 


-8 ± 5 


2 





Kaminski et al. 
TABLE 3 — Continued 



^obs 

(A) 


Alab 

(A) 


Molecule 


Electronic 
transition 




Feature 


Vel. 

(kms-i) 


Ref. 


5978.5 


5977.433 


TiO 


■y 


(3,0) 






3 


5980.8 




YO 




(0,0) 




-8 ± 5 


9 


oyoo.o 


oyo 1 . 1 z 


I w 


A2n , Y2y^+ 
^^3/2 ^ 






-i-f^si -t- ^ 


2 


oyyy.z 


OUUi.lol 


i IKJ 


7 


\o,\J ) 






Q 


6027.8 


6028 1 84 


TiO 


7 


(3,0) 


Ri 




3 


6034.4 


6036.17 


ScO 


A^n - X^E+ 


(0,0) 


^2Cr ^ ^*-2Cf 


-55 ± 10 


2 




6036.684 


TiO 


7 


(4*1) 






3 


6061.9 


6060.782 


TiO 


7 


(4,1) 


R2 


+55 ± 5 


5 


6065.2 


6064.31 


ScO 


A^n - x2s+ 


(0,0) 




+45 ± 10 


2 




6072.65 


ScO 


A^n - x2s+ 


(1,1) 






2 


6077.5 


6079.30 


ScO 


A^n - x2s+ 


(0,0) 


QQir- + "^Rir- 


-55 ± 10 


2 




6088.245 


TiO 


7 


(4,1) 


Ri 




3 


6116.0 


6115.97 


ScO 


A^n - x^s+ 


(111) 






2 


OIOU.O 


u±oz.uy 1 


YD 
I w 


A2tt -v2y^ + 

^-^1/2 ^ 




W21 




Q 


6135.0 


6134.472 


TiO 


7' 


(0,0) 


To 
K3I 




3 


6144.5 




TiO 


7' 


(0,0) 








6146.9 


6148.68 


TiO 


7' 


(0,0) 


■5R21 


-90 ± 10 


6 


6157.1 


6158.52 


TiO 


7' 


(0,0) 


Ri 




6 


6185.1 


6186.32 


TiO 


7' 


(0,0) 


R2 




6 




6189.65 


TiO 


7' 


(0,0) 


Q2 




6 


6213.4 


6214.93 


TiO 


7' 


(0,0) 


R3 




6 


6220.9 


6222.26 


TiO 


7' 


(0,0) 


Q3 


-70 ± 10 


6 


6269.2 


6268.86 


TiO 


7' 


(1,1) 


R3 


± 15 


5 


6276.3 


6276.03 


TiO 


7' 


(1.1) 


Q3 


± 15 


5 


6293.2 


6294. 8U 


TiO 


7 


(2,0) 


R3 




5 


6319.4 


Do2i.2i 


TiO 


7 


(2,0) 


R2 


-77 ± 5 


5 


6349.5 


Q c; 1 on 


1 lU 


7 


(2,0) 


Hi 


—80 + 5 


5 




DOO/ .OO 


1 lU 


7 


(3,1) 


H3 




5 


6384.5 


6384.18 


TiO 


7 


(3,1) 


R2 


+58 ± 5 


5 


6410.2 


6408.26 


ScO 


A 9T-r -^7-9^4- 

A^n3/2 - x^s+ 


(0,1) 


"Q2G + "R2G 




2 


6415.1 


6414.75 


TiO 


7 


(3,1) 


Rl 


+58 ± 5 


5 


6459.6 


6457.78 


ScO 


A2ni/2 - X2S + 


(0:1) 


'^QlG + '^RlG 




2 


6448.9 




TiO 


7 


(4,2) 


R2 






6649.3 


6651.271 


TiO 


7 


(1-0) 


R3 




1 


6679.1 


6680.796 


TiO 


7 


(1,0) 


R2 




1 


6712.5 


6714.477 


TiO 


7 


(1,0) 


Rl 




1 




6717.599 


TiO 


7 


(2,1) 


R3 




1 


6747.0 


d747.d13 


TiO 


7 


(2,1) 


R2 




1 


Dr<5i.i 


D ( cSi.cSiO 


1 IKJ 


7 


(2,1) 


Kl 




1 






TiO 


7 




rt3 




\ 


6815.9 


6815.139 


TiO 


7 


(3 2) 


R2 




1 


6832.2 


6832.596 


TiO 


bin - X^A 




R 


-82 ± 8 


3 


6838.9 


6838 


TiO 


b^n - X^A 




Q,P 




3 


6849.9 


6849.927 


TiO 


7 


(3,2) 


Rl 




1 


6884.3 


6885.25 


VO 


B^n - x*s- 


^2 0^ 


^R42 




2 


6SQ9 n 


uoyo.oyj 


VO 


B'*n — X'*s~ 




^42 




2 


oyoo.o 


Dyo ( . 


V w 


r4tt "Y"4v— 
11 — .A. Zj 








6950.9 


6952.39 


VO 


B ll — X 2j 


(2,0) 








6977.9 


6978.51 


VO 


B"*!! - X*S~ 


(2,0) 








6986.6 


6985.59 


VO 


B*n - x*E- 


(2,0) 








7003.5 


7003.40 


TiO 


7 


(0,0) 


■'R32 




4 


7037.8 


7038.28 


TiO 


7 


(0,0) 


■5R21 




4 




7054.256 


TiO 


7 


(0,0) 


R3 




1 


7085.5 


7087.566 


TiO 


7 


(0,0) 


R2 




1 




'71 OA norv 


liU 


7 


(1,1) 


R3 




1 


/ izo.4 


71 oc^ c^i n 

( 1.^10. OlU 


1 IKJ 


7 




rti 




1 






TiO 


7 


(1 l\ 


Ro 




1 


7231.8 


7230.774 


TiO 


7 


(2,2) 


R2 




3 


7270.9 


7269.985 


TiO 


7 


(2,2) 


Rl 




3 


7332.5 


7333.93 


VO 


B^n - x^E- 


(1,0) 


'^R42 


-60 ± 8 


7 


7371.5 


7372.70 


VO 


B*n - x*s- 


(1,0) 


■^Qsi 


-60 ± 8 


7 




7381.70 


VO 


B^n - x^s- 


(1,0) 


«P31 




7 




7393.34 


VO 


B^n - x^E- 


(1,0) 


^R21 




7 




7405 


VO 


B^n - x^E- 


(1,0) 






7 




7416.13 


VO 


B^n - x^E- 


(1,0) 


'3P21 




7 




7433.79 


VO 


B^n - x^E- 


(1,0) 


Rl 




7 




7442.9 


VO 


B^n - x^E- 


(1,0) 






7 


7453.8 


7455.2 


VO 


B^n - x*E- 


(1,0) 






7 




7474 


VO 


B^n - x*E- 


(2,1) 




+58 ± 8 


7 


7589.8 


7589.284 


TiO 


7 


(0,1) 


R3 




1 



g 
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-^obs 

(A) 


-^lab 

(A) 


Molecule 


Electronic 
transition 




Feature 


Vel. 

(km s"-"- 


) 


Ref. 


Notes 




'TKOT TOO 


i IKJ 


7 


1^1 1 J 








1 
1 




7704.2 


77nzL QQ7 


TiO 


7 


\^^^} 


-n.2 






1 

-L 




7744.4 


7743.056 


TiO 


7 


(2,3) 


Ra 






1 




7750.9 


7749.487 


TiO 


7 


(1,2) 


Ri 






1 




7783.8 


7782.854 


TiO 


7 


(2,3) 


R2 






1 




7863.20 


7865.29 


VO 


B^n - x^s- 


(0,0) 


■5Q42 


-80 ± 


4 


8 


h 


7894.20 


7896.38 


VO 


B^n - x^s- 


(0,0) 


^Qai 


-80 ± 


4 


8 


h 


7905.72 


7907.80 


VO 


B^n - x-^S' 


(0,0) 


«P31 


-80 ± 


4 


8 


h 




7919.12 


VO 


B^n - x^E- 


(0,0) 


^R21 






8 




7937.55 


7939.62 


VO 


B-^n - x^E- 


(0,0) 


R2 


-80 ± 


4 


8 


h 




7950.29 


VO 


B^n - x^s- 


(0,0) 


'5R23 






8 





References. — fl) [Phllilpa ifTOTt); (2) IKopp et all (ITotI): (3) ISchwenkd (11991): (41 |Ple3 ifTOOl): (5 ) IJflrgensenl |T99i ): 
(6) IHocking et al] (|1979l '): (Tl lCheung et all (|1994l '): (Sl IAdam et all l|1995l ): (9) IBernard fc Gravinal (|198a ). 

After correcting the laboratory wavelength from lPhillipsI (|1973f ) for the isotopic shift, as found from data in ISchwenkd l|1998f ). 
one gets for the isotopomer *®TiO the predicted stationary wavelength of 4759.205 A and the velocity of -80kms~^. 
^ The bandhead is strongly contaminated by the rotational structure of the a (2,0) band of TiO in the range 4793-4803 A. 

The bandhead is blended with two Fe I absorption lines at 5166.286 A and 5168.901 A. However, the blue edge of the blend 
should be formed by molecular absorption. The velocity given in the table is for the *^TiO isotopomer (laboratory wavelength 
concluded from ISchwenkdfl998l ) . 

All databases give different laboratory wavelengths. 
^ The band is probably contaminated by the A'^ni/2-X^E'*' (1,0) Qi feature of ScO. 
^ In emission. 

s Velocity determined in simulations, see Fig. [14] for details. 

^ In the table the velocity of the blueshifted component of VO is given. The fit to all features of the VO (1,0) and (0,0) bands 
indicates a velocity of -2kms~^. 



